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Separation of instrumental and interstellar components 

CP. de Vries^, J.W. den Herder^, J.S. Kaastra^, F.B. Paerels^, A.J. den Boggende^, and A. P. Rasmussen^ 

^ SRON, National Institute for Space Research, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands 

^ Columbia Astrophysics Laboratory, Columbia University, 550 West 120"" Street, New York, NY 10027, USA 

Received -< datel >-, accepted -< date2 >~ 

Abstract. High resolution X-ray spectra of the Reflection Grating Spectrometer (RGS) on board the XMM satellite 
are used to resolve the oxygen K absorption edge. By combining spectra of low and high extinction sources, the 
observed absorption edge can be split in the true interstellar (ISM) extinction and the instrumental absorption. 
The detailed ISM edge structure closely follows the edge structure of neutral oxygen as derived by theoretical 
R-matrix calculations. However, the position of the theoretical edge requires a wavelength shift. In addition the 
detailed instrumental RGS absorption edge structure is presented. All results are verified by comparing to a subset 
of Chandra LETG-HRC observations. 
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1. Introduction 

In the past, detectors in the low energy (0.3 - 2.0 keV) 
X-ray band offered spectra of modest spectral resolution. 
To properly model these spectra, equivalent spectral res- 
olution knowledge of the oxygen K-edge absorption was 
required in order to correct the intrinsic source spectra 
for the effects of interstellar extinction. With the appear- 
ance of high resolution X-ray spectrographs, notably the 
XMM-Newton RGS (den Herder et al. 120011) and Chandra 
LETG (Brinkman et aL l^nnUjl . the ISM oxygen K-edge ab- 
sorption can be studied in much greater detail. These in- 
struments, however, usually suffer from absorption by oxy- 
gen located within the instruments. By comparing spectra 
from sources with low and high extinction and assuming 
that the interstellar cross-section per hydrogen atom is 
constant over the sky, the instrumental and interstellar 
components can be separated. 

Models of extinction in X-rays by the interstel- 
lar medium have been constructed, amongst others, by 
Morrison and McCammon H1983fl . These models are based 
on assumed abundance's of the ISM and cross sections 
by Henke et al. (1982). However, the spectral resolution 
in current spectrographs surpasses the resolution of the 
Henke tables. 
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For neutral atomic oxygen far more detailed absorp- 
tion cross sections have been calculated by the R-matrix 
approach (McLaughlin and Kir by, 11998(1 . These cross sec- 
tions include all the fine structure resonance's along the 
K edge (around 22.6 - 22.9 A), including the ls-2p transi- 
tion at approximately 23.5 A. Neutral oxygen is thought 
to be the major constituent of the ISM. Therefore, a good 
model of the detailed spectral structure around the oxygen 
K edge can be constructed by combining the "R-matrix" 
cross sections of neutral oxygen with the Morrison and 
McCammon extinction cross sections for the remaining 
elements. The ISM oxygen abundance as presented by 
Wilms et al. |2QQQ|) is used. Fig. ^ shows the combined 
cross section. 

In the matrix calculations uncertainties of up to 3.5 eV, 
or 0.15 A are reported for certain energy levels which can 
lead to shifted absolute positions of the resonance struc- 
tures. The relative structure however, which gives rise to 
the detailed shape of the observed edge, can be expected 
to be known with much higher accuracy. The theoretical 
bound-bound resonance's wavelengths (e.g. the ls-2p line) 
were shifted to be in agreement with laboratory data, such 
that for these transitions the error can be expected to be 
less. 

Of course interstellar oxygen can also be present in 
other forms then neutral atoms. Ionised oxygen can be 
expected and oxygen can be bound in molecules (e.g. OH, 
CO, H2O), either as free gases, or as a solid component 
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Fig. 1. Theoretical cross sections for the ISM around the 
oxygen K edge based on the R-matrix (McLaughhn and 
Kirbv. [T998|l computations for neutral atomic oxygen com- 
bined with the Morrison and McCammon (jliJS^ profile 
for other elements. The edge is not infinitely sharp, but 
consists of a series of resonance features, indicated in the 
plot. The well known ls-2p transition feature is also indi- 
cated. 

of interstellar dust (e.g. ice, sihcates). Although the X- 
ray signature of ionised oxygen is well known, the X-ray 
signature of bound oxygen is much harder to identify. 
Depending on the detailed chemical and crystal structure 
and the presence and structure of other compounds, the 
X-ray spectral signatures may shift by up to 5 % (chem- 
ical shift, see e.g. Sevier, I1979|l . and can be broadened. 
For this reason the exact X-ray absorption spectrum of 
interstellar dust is not really known. Since neutral atomic 
oxygen is believed to be the most significant constituent 
of interstellar oxygen it is therefore justified to include its 
X-ray absorption features only in our interstellar oxygen 
model. 

In this paper we derive the shape of the interstellar 
oxygen absorption edge by comparing the absorption fea- 
tures in sources with low and high extinction. We com- 
pare the neutral oxygen model cross sections with the de- 
rived shape of the spectral profiles of interstellar extinc- 
tion observed by the RGS, corrected for the instrumental 
response. Possible deviations from the model will be dis- 
cussed in terms of the presence or upper limits of other 
oxygen components. In addition the detailed features of 
the RGS instrumental absorption are resolved. A subset 
of Chandra measurements is used to crosscheck the results 
obtained. 

2. Data 

To study spectral absorption features of the oxygen K 
edge, a number of bright sources with simple intrinsic 
(power law) spectra and no obvious narrow spectral fea- 
tures, were selected. Table ^ lists the sources. The cor- 
responding XMM-Newton orbits of the RGS data and 



Chandra observation ID's of the LETG data are shown, 
as well as the A^h ISM column densities used. 

Since the main features in our spectra are caused by 
oxygen, we also show the oxygen column density based 
on the ISM abundance by Wilms et al. H2()()();i . For the 
extragalactic sources (Mrk 421, PKS 2155-304, 3C 273) we 
checked that we got reasonable agreement with the oxygen 
spectral features, starting from the Nyl columns from an 
independent source fKaastra. l2()()2|l . while for the galactic 
sources (Sco-Xl and 4U0614-I-091) we actually fitted the 
oxygen column and convert to A^h- It may be noted that 
for the derivation of the actual shape of the interstellar 
oxygen edge the precise oxygen column densities of the 
individual sources are irrelevant. 

Finally, for comparison only, we add the A^h column 
densities based on Hi radio surveys. Since the line of 
sight towards the galactic sources is likely to cross galac- 
tic molecular material, a correction for the presence of H2 
based on a ratio of H2 to H i of 0.3 (Dame et al. ll987|l is ap- 
plied. For Sco-Xl we assume all observed H l to be in front 
of the source, since its distance of 2.8 kpc (Bradshaw et al. 
fro99f) and galactic latitude of about 24 degrees places it 
well above the scale height of « 150 pc of H i in the galaxy. 
For 4U 0614-1-091 a wide variety of distances is reported 
(see Brandt et al.. ll992|l from 2 to 8 kpc although an up- 
per limit of 3 kpc seems likely. Combined with the low 
galactic latitude for this source of only 3 degrees, makes 
the estimated Hi column highly uncertain. 

Despite some uncertainties in the Hi column densi- 
ties, tabled shows quite good agreement between our A'h 
columns derived from X-ray data and the Hi columns. 
This is a strong indication that the absorption is truly in- 
terstellar and not originating at the source. This means 
the absorption characteristics can be assumed similar for 
all our sources. 

On the RGS the oxygen edge is covered by CCD de- 
tector 4. Unfortunately, early in the mission the CCD4 
chain on RGS2 failed. As a consequence the oxygen data 
for on-axis sources originate from the single CCD chip on 
RGSl only. The bad spots (warm or hot pixels) on this 
CCD detector 4 cause large systematic uncertainties in the 
spectra so these regions have to be omitted from the analy- 
sis. These bad spots usually occur at the same wavelength 
for one source, but could be slightly shifted for a different 
source, due to small differences in source position with re- 
spect to the telescope axis. By combining the spectra of 
the low extinction sources Mrk 421 and PKS 2155-304 the 
efi'cct of these bad spots is reduced significantly due to 
slight differences in the relative pointings of these sources. 

For the "high extinction" source Sco-Xl a different ap- 
proach was adopted. Sco-Xl is sufficiently bright to allow 
large offset pointings (of about 27 arcmin.) which shifted 
the oxygen edge to CCD detector 3 (orbit 402). Since both 
CCD3 detectors are active on both RGS's, redundancy 
between the two CCD3 detectors solves the problem with 
bad positions on the detector corresponding to a particular 
wavelength range. Comparing the individual CCD3 spec- 
tra and an additional on-axis pointing on Sco-Xl with the 
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Table 1. Overview of sources and basic data 





Observation 


exposure' 


A''h column 


No column^ 


A''h from H I 


Source 


XMM-Newton ^ 


Chandra 


(ksec) 


(xlO^^cm-^) 


(xlO'^cm-2) 


(xlO^^cm-^) 


Mrk 421 


84,171,259 




172 


1.66 ±0.09"' 


0.85 ±0.04 


1.61 ±0.10'= 


PKS 2155-304 


87,174,362 


3166 


346 


1.25 + 0.08"^ 


0.65 ±0.04 


1.36 ±0.10" 


3C 273 


277 




90 


1.76 + 0.10^^ 


0.90 ±0.05 


1.68 ±0.10"= 


Sco-Xl 


224,402 




12 


20.0 ± 1.5^ 


10.2 ±0.8 


19.0 ± 3^ 


4U 0614+091 


231 


100 


23 


31. ±4.'' 


15.9 ±0.4 


28 ± 10*5 



Derived from the A^h column depth based on the oxygen ISM abundance's from Wilms, Allen and McCray ()2()()()^ 
^ XMM-Newton orbit number. 
" Chandra Obs-ID. 

Taken from independent Chandra data (Kaastra, 120021 
" Based on Hi observations by Lockman and Savage Ill995t . 
' Fitted on the RGS data (this paper). 

^ Based on the Hi survey by Hartmann and Burton (119971 and a ratio of H2 to Hi of 0.3 (Dame et al. 119871 . The errors are 

based on local gradients in the H I map and uncertainties in source distances. 
^ Based on fits in both the Chandra and RGS data (this paper). 
' Total RGS exposure time 



oxygen edge on CCD4 (orbit 224) it was confirmed that 
the instrumental oxygen edge is identical for the different 
detector chips, which implies that the same instrumental 
profile applies irrespective of the location of the Oxygen 
edge on the instrument. 

The data of 3C 273 were used as an independent check 
to verify that the derived instrumental absorption profile 
for the RGS is correct. 

In all figures in this paper data are fluxed spectra, 
corrected for exposure and the effective area excluding of 
course the effects of instrumental oxygen absorption. Error 
bars represent the la statistical error. 

The lower effective area and limited exposure times 
of Chandra imply higher statistical noise compared with 
XMM. Including the Chandra observations however, en- 
ables an independent check of the derived properties of the 
interstellar oxygen absorption. All Chandra data shown 
are corrected for the overlapping 2"^^ and orders. 

3. RGS analysis 

Fig.[21shows the RGS spectrum around the oxygen K edge 
for Sco-Xl. Overplotted is the theoretical profile based on 
the cross sections of Fig.^ convolved with the RGS instru- 
ment response and using an extinction column depth of 
Nn = 2.0 X 10^1 cm~2. With this column depth the ls-2p 
absorption line fits the data after applying a 30 mA shift 
of the theoretical position to longer wavelengths. This is 
somewhat larger than expected based on the known ac- 
curacy of the RGS wavelength scale of Icr = 8 niA (den 
Herder et al.. \Mn^ . 

Due to the very high optical depth of the narrow Is- 
2p line core, saturation effects can be expected. Given the 
cross sections of the ls-2p line as shown in Fig. ^ the ex- 
pected relation between equivalent width of the line and 
the column density was computed for a variety of veloc- 
ity dispersions. This is shown in Fig. |3 The bend in the 
curves for high column density does indicate saturation 
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Fig. 2. RGS spectrum around the oxygen K edge of Sco- 
Xl (grey data points). The solid black line represents the 
theoretical interstellar absorption (Fig. ^ convolved with 
the RGS instrumental resolution for an ISM column den- 
sity of Nu = 2.0 X lO^^cm^^^ ^^j^g jgj^ 2s-2p absorption 
feature is also indicated. 

effects, but due to the broad wings in the line (contrary 
to e.g. lines of O i in the UV), the equivalent width of the 
line keeps increasing steadily even for high optical depths. 
This means that even for high extinction the column den- 
sity can be reasonably well determined from the equivalent 
width. However, the actual equivalent width and depth of 
the observed line profile do depend on the velocity disper- 
sion along the line of sight. Entering realistic velocity dis- 
persions (< 200 km/s) did not improve the fit sufficiently 
to allow for tight constraints on the velocity dispersion and 
column depth. Therefore a fit of the ls-2p line only, leaves 
a rather large uncertainty on the actual column density. 
Fixing the column depth to A^h = 2.0 x 10^^ cm^^ and 
performing more detailed modeling of the ls-2p line depth 
using the SPEX package (Kaastra. I2()()2c|l points to a ve- 
locity dispersion in the range of 30 - 50 km/s. However, 
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Fig. 3. Relation between equivalent line width of the ls-2p 
absorption and A^o column density. The solid lines show 
the curves for neutral oxygen for different velocity dis- 
persions of 10, 30 and 100 km/s from the lowest to the 
highest solid line respectively. Curves for velocities < 10 
km/s overlap with the lowest curve (10 km/s). The broken 
lines show the same for singly ionised oxygen (Oil). 



this is highly dependent on the real distribution of the 
velocity components. In addition, given the errors on the 
equivalent widths (TableEJ and too poor resolution of the 
RGS for this purpose, the fitted velocity dispersion can 
only be seen as an indication only. 

Since the depth of the absorption below the edge (be- 
low 22.5 A) is a measure for the total of ISM and instru- 
mental extinction, it was checked whether the difference 
between ISM absorption model and data below the edge 
can be attributed to the instrument. The depth of the in- 
strumental oxygen absorption determined later in this pa- 
per, does indeed match this offset. Since the depth of the 
absorption below the edge extends over a large wavelength 
range and is therefore not influenced by ISM velocity dis- 
persions, it is believed that the current fit does give a good 
measure of the true oxygen column depth for this source. 
The derived column density is also in good agreement with 
the numbers retrieved from Hi observations (see Table 

Fig.^shows the RGS spectrum for the sum of Mrk 421 
and PKS 2155-304 data. This is the "low extinction" data. 
Of course this summed spectrum yields a kind of weighted 
average of the spectral characteristics of the two individ- 
ual objects Mrk 421 and PKS 2155-304. Since we are not 
interested in the intrinsic source spectra the large scale 
shape of the combined spectrum is not relevant. The A^h 
column depth for both sources is low and does not differ 
too much between the sources. It is therefore estimated 
that the distortion of the ISM absorption profile due to 
adding the two sources is well within the statistical uncer- 
tainty of the individual data points. 

For these "low extinction" data, the model wave- 
lengths are shifted by 35 mA to longer wavelengths to 
fit the ls-2p feature. This shift is consistent with the 
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Fig. 4. RGS combined spectrum around the oxygen K 
edge for the low extinction sources (Mrk 421, PKS 2155- 
304). A fit of the theoretical ISM absorption model to 
the ls-2p line at 23.5 A is also shown (black line). Main 
difference with the data points is due to the additional in- 
strumental absorption. Two key instrumental features at 
23.05 and 23.35 A are indicated. 



30 mA shift in Fig. [21 given the 8 mA error on the RGS 
wavelength scale. Main difference between the data points 
of these low extinction sources and the overplotted ISM 
absorption is due to the additional instrumental absorp- 
tion. Two key instrumental features can be clearly recog- 
nised at around 23.05 and 23.35 A. This is the fingerprint 
of the actual compounds in the instrument which cause 
the absorption. These dips can also be noticed in Fig. El 
in the Sco-Xl data. The fact that the equivalent width of 
the 23.35 feature (see Tablej^J is the same for both Sco-Xl 
and the low extinction sources, despite the fact that the 
A^H column densities differ more than a factor of 10, does 
confirm the instrumental nature of this feature. 

Taking the ratio between the "low" and "high" extinc- 
tion data will remove all the features which are constant 
between the data, in particular the instrumental features. 
What is left must represent the true interstellar extinction. 
This ratio is shown in Fig. (black data points). 

The theoretical ISM extinction curve (broken grey 
line) is fitted to the depth of the edge below 22.5 A. By 
comparing the theoretical oxygen edge with the extinction 
data, it appears that the shape of the edge matches the 
data, but that the position of the edge is shifted. The solid 
grey line shows this effect. It is identical to the broken line 
but with the edge (all the resonance structures identified 
with the indicated edge in Fig. shifted by 0.13 A, to 
larger wavelengths. The error uncertainty of this shift is 
estimated at 0.02 A. Hereafter, we use these absorption 
cross sections with shifted oxygen edge as our model of 
the interstellar oxygen K edge. 

It shows clearly in Fig. that the absorption dips 
around 23.05 and 23.35 A, which are present in both 
the low and high extinction data (Figs. El and have 
largely disappeared in the ratio plot, which indicates that 
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Fig. 5. Ratio of the spectra of Figs. |21 and 0] (black data 
points). This represents the variable interstellar extinction 
between the sources. The broken grey line shows the fitted 
theoretical extinction. The solid grey line shows this the- 
oretical extinction with the edge shifted by 0.13 A. Error 
bars shown represent the la statistical error. 



both features are generated by the RGS instrument itself. 
Possible residues deviate less than 3ct from the model. 

Any deviations in Fig. [S] from the model can point 
to absorption features of other components, in particu- 
lar other forms of oxygen. Any amount of oxygen, either 
atomic, in molecules or solid compounds will be reflected 
by the depth of the oxygen edge. Therefore the depth of 
the edge by itself cannot separate different components. 
Oxygen compounds other then the neutral atomic form 
can only be seen if they show distinct, sufficiently narrow, 
absorption features. In order to judge possible strengths 
of other absorption features, the level of the fitted contin- 
uum is important. In Fig. [S] it shows that the continuum 
is very well constrained by the data below 22.5 A and 
above 23.7 A. Although fit residues between 22.5 and 
23.5 A seem to suggest additional structure, the tight con- 
straints on the continuum level do not allow much free- 
dom to move the continuum up to a higher level such that 
there is room for clear absorptions. Taking into account 
the error on the fitted continuum, the fit residues can be 
analysed in terms of upper limits on possible other oxygen 
compounds. 

A striking feature in the fit residue is the feature 
around 22.77 A. At this wavelength some excess fiux ap- 
pears in the ratio plot (Fig.[SJ, which means that this fea- 
ture is not present in the Sco-Xl data (Fig.|51). Inspection 
of the individual XMM data show that this feature is 
present in both Mrk 421 and PKS 2155-304 data but not 
in the 4U 0614-f91 and Sco-Xl spectra. The spectrum 
of the combined low extinction sources. Fig. shows this 
feature as a clear distinct absorption feature. It is unlikely 
that the apparent excess fiux in the ratio plot is caused by 
two absorption dips on both sides of the 22.77 A feature, 
e.g. due to solid ISM oxygen compounds. In that case the 



Sco-Xl should have shown the 22.77 A feature in "emis- 
sion" as well. In the subsequent section we will show that 
the Chandra data of PKS 2155-304 and 4U 0614-^91 show 
the same behaviour as the RGS data. For this reason this 
feature is attributed to the PKS 2155-304 and Mrk 421 
source spectra. 

If we tentatively identify this feature (Fig. ^ in the 
PKS 2155-304 and Mrk 421 spectra with an Oiv ab- 
sorption blend at 22.75 A, given the equivalent width 
of « 10mA, we derive a column density of Nqw of 6 x 
lO^^cm^^. The most likely candidate for this absorption 
is the local filament as observed e.g. in Ovil by Nicastro 
et. al. l|^02|i) . with a column density of 4 x 10 ^^cm"^. 
However, since Nicastro et. al. derive an Ovi column of 
only 0.1 X lO'^^cm"^, the Oiv is most likely due to a less 
ionised component. If the gas is photo ionised, we expect 
Oiv to peak at log^ ~ —0.8. For this ionization parame- 
ter the column density of the neighbouring O ill and O v 
ions is a factor of 3 smaller than the O IV column and we 
estimate the strongest spectral features of those ions to be 
at 22.33 A(Ov, eq. width = 4 mA) and 23.08 A(Oiii, eq. 
width = 2 mA) . Equivalent widths of this order of mag- 
nitude are not completely excluded by the data, taking 
into account that these lines positions are uncertain by 
about « 0.03 A. Unfortunately, there is no clear positive 
indication of these lines. 

Looking in detail to other features in the fit residues 
of the ratio spectrum of Fig. |S1 features can be recognized 
at 23.05 and 23.35 A. The 23.05 A feature seems an in- 
complete cancellation of the 23.05 A instrumental feature 
and has an equivalent width of 3 ± 3 mA. Finally there 
seems some structure around 23.35 A which has an equiv- 
alent width of 4 ± 3 mA. A tentative explanation could be 
traces of O iii (at 23.05 A) and O ii (at 23.35 A). This im- 
plies column densities of « 5 x lO^^cm^^ (see Fig.OJ. All 
these column densities are sufficiently low, such that their 
combined K-edge has an optical depth less than 0.01. This 
puts constraints on the average ionization stage of oxy- 
gen in the interstellar medium. Higher ionization stages 
of Oxygen may dominate the hotter components of the 
interstellar medium. 

Apart from possible lines of ionised oxygen, oxygen 
locked in solids may cause features in the fit residues of 
Fig. 1^1 It is well known that in solids the spectral features 
of X-ray absorption lines can be shifted (see e.g. Sevier, 
I1979|l and broadened with respect to the corresponding 
lines in free atoms. For broader lines saturation effects 
are less critical. If the combined features identified above 
would be attributed to shifted and broadened ls-2p lines 
of neutral oxygen in solid compounds, we estimate the 
amount of oxygen in these features at less than 10 % of the 
neutral oxygen. In case of possible saturation effects the 
estimated upper limits will become smaller. If however the 
combined chemical shifts of lines of oxygen locked in solids 
smear out its signature over a broad wavelength band, 
its signature will be completely lost and large amounts 
of oxygen may be hidden in solids which are simply not 
detected by X-ray absorption features. 
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Fig. 6. Fit of modified ls-2p lines to the instrumental fea- 
tures for the low extinction data (Fig.0}. The broken line 
shows the combination of our interstellar extinction model 
(black line) and the fit of a broadened shifted atomic ls-2p 
line model to the 23.05 and 23.35 A instrumental features. 
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Fig. 7. The derived instrumental efficiency due to the in- 
strumental oxygen absorption. The extrapolations beyond 
the wavelength range of 21.8 - 24.0 A have been derived 
from the Henke H1982|l tables (broken lines). 



The fact that the depth of the ls-2p line does not quite 
match the data is not really relevant. As explained before, 
the relative depths of the edge and the line of the high 
extinction component in the ratio data is rather arbitrary 
and depends on other factors (ISM velocity field) then just 
the column depth of the absorbing material. 

If the instrumental features at 23.05 and 23.35 A are 
interpreted as shifted ls-2p lines of bound oxygen in 
some solid compounds located on the detectors or op- 
tics, the column density for these features can be fitted. 
First our interstellar extinction model is fitted to the in- 
terstellar ls-2p line of the low extinction data at 23.5 
A(Fig. EI). Fig. El shows that for the fitted ISM column 
density of A^h = 1-5 x 10^" cm^^, which corresponds 
to No = 0.76 X 10^^ cm~^, line saturation effects are 
minor and only uncertainty in velocity dispersion will in- 
fluence the determined column density. We assumed veloc- 
ity dispersions < 10 km/s for these high galactic latitude 
sources. 

In addition to this ISM absorption the instrumental 
23.05 and 23.35 A features are fitted by shifting and broad- 
ening the ls-2p line. It now appears that with this instru- 
ment model the depth of the edge below 22.5 A is also well 
fitted. This is interpreted as a confirmation that the 23.05 
and 23.35 A features are indeed shifted and broadened 
ls-2p lines. The oxygen column density in both instru- 
mental features amounts to A^o — 2.5 ± 0.2 x 10^^ cm~^. 
If all this instrumental oxygen would be in the form of 
water(ice) only, this column density amounts to a layer of 
about 750 ± 50 A thickness. However, the fact that we see 
two features might indicate that the oxygen is present in 
two different compounds. 

The actual difference in Fig. between the fitted theo- 
retical interstellar extinction curve and the data (ignoring 
the 22.77 A feature) can be interpreted as the instrumen- 
tal oxygen absorption. For these low extinction sources, 
any uncertainty in the exact shape of the interstellar oxy- 



gen edge, whether or not features of other forms of oxygen 
then just neutral oxygen are present, hardly plays a role. 
The derived corresponding efficiency of the RGS instru- 
ment is plotted in Fig. [7| This efficiency is implemented in 
the current calibration files (CCF) for the RGS response, 
used by the XMM SAS data reduction software package 
(Jansen et al. .2001^ . Note that some of the sharp "edges" 
in the graph result from the fitting process used to derive 
the instrumental curves. They have no physical meaning. 
These sharp features disappear when this correction is ap- 
plied to the instrument response. 

Our interstellar and instrumental absorption models 
around the oxygen K edge were verified by applying 
them to 3C 273. The statistical and systematic (due to 
warm/hot columns in the RGS CCDs) noise is higher for 
this data than the datasets used above. Using a column 
depth of Ah = 1.76 x 10^° cm~^, derived from indepen- 
dent Chandra data, yields an acceptable fit to these data 
(Fig.©. 

4. Cross checking on Chandra data 

The method to compare data from low and high extinction 
sources to derive absorption characteristics of the ISM can 
also be applied to Chandra data. We looked at a subset 
of Chandra LETG HRC data, since this Chandra grat- 
ing/detector combination has the highest effective area at 
the oxygen K-edge. The observations are listed in Tabled 
Fig. El shows the Chandra observation of the low ex- 
tinction source PKS 2155-304. Comparing with the XMM- 
RGS observation of the low extinction sources in Fig. 01 it 
is immediately apparent that the deep instrumental ab- 
sorption feature seen at 23.05 A in the XMM data is 
largely absent in the Chandra data. The 23.35 A instru- 
mental absorption feature however can be recognised in 
the Chandra data. This feature is treated similarly as for 
the XMM-RGS, by assuming it to be a shifted/broadened 
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Fig. 8. Fit of our interstellar model with a column depth 
of Nn = 1.76 X 10^°cm~^ and our instrumental efficiency 
model to the RGS spectrum of 3C 273 in the wavelength 
range around the oxygen K edge 
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Fig. 9. Chandra LETG HRC observation of PKS 2155- 
304. The solid black line represents the ISM absorption 
using our iVn of 1.25 x lO^^cm"^. The broken line shows 
a fit to the instrumental line at 23.35 A which can also 
be recognised in the XMM-RGS data. Here, this fit corre- 
sponds to an oxygen column depth of 8 x lO^^cm"^. The 
feature at 22.77 A, which is also seen in the XMM-RGS is 
attributed to the source itself. 
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Fig. 10. Chandra LETG HRC observation of 
4U 0614-1-091. The sohd line shows the modeled in- 
terstellar absorption using the Chandra wavelength 
resolution, while the broken line adds the instrumen- 
tal absorption as determined in the PKS 2155-304 
observation (Fig. El- 
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Fig. 11. XMM-RGS observation of 4U 0614+91. These 
data can be directly compared with Fig. ^] The broken 
line adds the RGS instrumental absorption. It can be seen 
that the ls-2p line in the Chandra data shows deeper ab- 
sorption due to the higher spectral resolution of Chandra. 



instrumental ls-2p line. This fit is shown with the plot- 
ted broken line. Also the suspected absorption due to the 
source at 22.77 A can be recognised in the Chandra data. 

From these observations it follows that the general in- 
strumental absorption characteristics between Chandra- 
LETG-HRC and XMM-RGS differ, but the 23.35 A com- 
ponent is common in both instruments. 

A high extinction source observed by Chandra is 
4U 0614-f091 (Paerels et al., WM^. The Chandra spec- 
trum of this source is shown in Fig. 1101 The correspond- 
ing XMM-RGS spectrum is shown in Fig. ^2 It can be 
noted that the Chandra ls-2p line shows a deeper ab- 
sorption than the XMM-RGS spectrum. This is due to 
the higher spectral resolution of Chandra of AA = 0.045 
FWHM (Chandra Proposers Observatory Guide) com- 
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Fig. 12. Ratio of the spectra of Figs-ITUlandl^fblack data 
points) . This shows the interstellar extinction according to 
Chandra. Overplotted (grey line) is our interstellar extinc- 
tion model which includes the 0.13 A shift of the oxygen 
edge. 

pared to AA = 0.065 FWHM (den Herder et al. [M7T|l 

for the RGS and means that this line is not resolved by 
the XMM-RGS. 

Both PKS 2155-304 and 4U 0614-1-091 Chandra data 
show an apparent wavelength inconsistency of 50 mA to 
larger wavelengths compared to the theoretical position 
of the ls-2p line. This can be compared with the 30 to 
35 mA shift in the same direction for the XMM-RGS 
data noted before. The systematic wavelength error on 
the Chandra data is estimated at 10 mA in our wavelength 
region (Kaastra et al.. l2002E)l . while we mentioned a wave- 
length inaccuracy of 8 mA for the RGS before. This means 
that, within the errors reported, the Chandra and RGS 
wavelength scales match, but that likely the R-matrix the- 
oretical ls-2p line position requires a slight adaptation 
of 35 ± 10 mA to longer wavelengths. This can be com- 
pared with the noted discrepancies between the theoreti- 
cal R-matrix computations and laboratory measurements 
of wO.25 % or 60 mA. Laboratory measurements report 
wavelengths between 23.49 A (Menzel et al.. I1996|l and 
23.53 A fStolte et aL. ITMTjl . 

In fig.^lthe broken black line shows the fitted model 
spectrum corrected for the chandra instrumental absorp- 
tion as derived from fig. |^ There is a hint for some ad- 
ditional (instrumental) absorption between 23.0 - 23.2 A, 
although the PKS 2155-304 data do not show this absorp- 
tion that clearly. However, the Chandra data are quite 
noisy in this respect. 

Combining the Chandra spectra of Figs. ^1 and 
1^ yields the interstellar extinction curve according to 
Chandra, shown in Fig. 1121 Since all instrumental fea- 
tures are removed from this plot, it can be directly com- 
pared with the corresponding XMM-RGS plot of Fig. [S] 
The solid grey line in both plots shows our ISM absorp- 
tion model which is identical in both Fig. [S] and Fig. [T^ 
Although the Chandra profile is much noisier, the pro- 



files are remarkably identical. Even the positive excursions 
around 22.77 A which hint at source absorption features, 
are repeated in the Chandra data although with far less 
significance. This similarity between the XMM-RGS and 
Chandra derived interstellar oxygen absorption features 
confirms the validity of the data analysis. 

5. Discussion and conclusions 

Comparing high resolution X-ray spectra of sources with 
low and high interstellar extinction, absorption features 
around the oxygen K edge can be separated into instru- 
mental and interstellar components. A summary of all fea- 
tures identified is presented in Table [51 

A clear signature of the RGS instrumental profile is 
seen through absorption structures at 23.05 and 23.35 A. 
These can be interpreted as the ls-2p line, broadened and 
shifted due to the bound state of oxygen on the instru- 
ment. A total column depth of oxygen on the RGS instru- 
ment of A'^o = 2 X 10^^ cm^^ is derived. 

The RGS instrumental oxygen edge is remarkably sta- 
ble over time. Early and later observations of Mrk 421 
(orbits 84 and 259) show, within statistical errors, exactly 
the same structure and magnitude of the edge. 

The 23.35 A feature is also seen in the Chandra LETG 
spectra. If the 23.05 and 23.35 A features point to differ- 
ent components, it is clear that the 23.35 A component is 
common in both instruments, while the 23.05 A compo- 
nent is less prominent in Chandra LETG/HRC. 

Likely compounds in the instruments which contain 
oxygen (apart from the known oxygen in the Chandra 
filter polyamide carrier) are water (ice) and oxidised metal 
surfaces. 

Both XMM-RGS and Chandra LETG have metal sur- 
faces, which will be slightly oxidised. One could therefore 
argue that the 23.35 A feature corresponds to metal ox- 
ides. For RGS, the active Si of the RGS CCD detectors is 
covered with a 270 A MgF2 insulation layer which, on top 
of that, holds an Al light shield, with a thickness in be- 
tween 450 and 740 A, depending on the CCD (Den Herder 
et al. l2()()l(l . In the Chandra HRC there is a separate filter 
which holds an aluminium layer. In both instruments the 
oxides can be resident on the aluminium filters or optical 
components. For RGS the computed oxygen column den- 
sity for the 23.35 A feature corresponds to an aluminium 
oxide layer of about 35 A. An independent measurement, 
using ellipsometry, of the oxide layer on the top of the 
Al light shield on one of the flight spare RGS CCDs in- 
dicated a layer of only w 10 A. Additional metal oxides 
will be present on the active Si as well on the optics. In 
the Chandra data the fit of the 23.35 A feature yields 
an equivalent aluminium oxide layer thickness of about 
110 A. 

The 23.05 A feature dominates the RGS instrument 
response. Contamination (ice forming on the cooled detec- 
tors) is very unlikely the origin of this absorption feature 
in view of the stable response of the instrument over a 
prolonged period of time during repeated observations of 
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Table 2. Summary of spectral features around the neutral oxygen-K edge 

wavelength Observed equivalent width^ (m^) in source spectrum 

(A) identification Mrk421 PKS 2155 3C273 Sco-Xl 4U 0614 

23.50 ±0.01'' neutral oxygen ls-2p inter- 16 ± 3 15 ± 3 20 ± 8 60 ± 5 70 ± 8 

stellar absorption line 20 ±5"* 80 ± 5"* 

22.6 - 23.0° interstellar oxygen-K edge )*■ )' Y 

absorption structure 

23.05 ±0.08"^ instrumental absorption )s )s )« )s )s 

23.35 ± 0.02 instrumental absorption 18 ± 4 15 ± 4 18 ± 6 19 ± 4 19 ± 3 

25 ±5'^ 30 ± 

22.77 ± 0.03 possible local intergalactic 9 ± 2 11 ± 6 <5 <3 <5 

absorption 

Equivalent width of feature in the RGS spectrum, unless noted otherwise. 
^ Shifted by 35 ± 10 mA to longer wavelengths with respect to the theoretical position of McLaughin and Kirby, ^1 99811 . 
Broad structure of Oxygen bound in a solid compound. The error also reflects the broad and asymmetric nature of the 
feature. 

Equivalent width of feature in the Chandra spectrum. 
" The K edge is not sharp, but extends over the indicated range. The edge is shifted by 0.13 ± 0.02 Awith respect to the 

theoretical R-matrix calculations. 
' An equivalent width has no meaning for an edge structure. 

^ Since this instrumental feature is very broad the continuum level and structure is very uncertain. For this reason no equivalent 
width was determined. 



Mrk 421. In addition, the relative intensities of the RGS 
aluminium and fluor calibration sources, which show X- 
ray lines at 1.487 and 0.677 keV respectively, did not 
change from before launch to recent times (orbit 440). 
These intensity ratios are very sensitive to accumulation 
of oxygen (water-ice) on the detectors. 

Another source of oxygen may be absorption of wa- 
ter by the MgF2 insulation layer during manufacturing of 
the CCD's, before the Al hght shield was applied. This 
process would lock the water behind the Al light shield, 
without any possibility to escape. MgF2 is known to be 
porous, and may be able to absorb considerable quantities 
of water, depending on the way it was deposited on the 
CCD surface. When all of the oxygen of the 23.05 A fea- 
ture would be in the form of water, we would need a layer 
of 550 A to be able to account for the observed column 
density. Despite the fact that the thickness of the MgF2 
layer is only known for the "solid equivalent" (assuming 
that the MgF2 is not porous), the required quantity of 
water seems quite large and is hard to account for solely 
in the 270 A MgF2. 

Laboratory measurements of RGS flight-spare CCDs 
are consistent with the instrumental oxygen edge mea- 
sured in flight. The systematic error for the laboratory 
measurements, however, is high compared to the flight 
measurement. 

The interstellar (ISM) extinction derived from XMM- 
RGS and Chandra LETG/HRC data is fully consistent 
between the instruments, confirming the validity of the 
analysis. 

The interstellar component is compared with theoreti- 
cal predictions based on the R-matrix (McLaughlin and 



Kirby, I1998|l computations of cross sections for atomic 
oxygen. There is indication that the theoretical ls-2p line 
position requires a slight shift of about 35 ± 10 mA to 
larger wavelengths. The structure around the K edge is 
shifted relative to the ls-2p line by 0.13 ±0.02 A to longer 
wavelengths (closer to the ls-2p line), when compared to 
the theory. The overall shape of the edge, which reflects 
the detailed resonance structure convolved with the in- 
strument response, fits the actual data well. The total 
interstellar X-ray absorption fits the expected profile for 
neutral atomic oxygen. No unambiguous additional signa- 
ture of oxygen bound in other compounds (e.g. interstel- 
lar dust grains) is found. However, we are only sensitive 
to sufficiently narrow features. If the type of solid state 
oxygen in dust grains would exhibit very broad absorp- 
tion features we would not recognise it. An upper limit 
of oxygen compounds other than neutral oxygen, which 
have distinct and narrow features around the oxygen edge 
is estimated at less than 10 % of the neutral oxygen. 

A possible interstellar absorption feature of the local 
intergalactic medium was found around 22.77 A in both 
PKS 2155-304 and Mrk 421 spectra. The exact compound 
causing this possible absorption feature is not unambigu- 
ous, but the fact that this feature is not seen in our galactic 
sources does point to an extragalactic origin. We tenta- 
tively identify this feature with O iv absorption. 
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